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Abstract 

With reference to the experimental observations by T. Yanagida and his co-workers 
on actomyosin interaction, a Brownian motor of fluctuating ratchet kind is designed 
with the aim to describe the interaction between a Myosin II head and a neighboring 
actin filament. Our motor combines the dynamics of the myosin head with a chemical 
external system related to the ATP cycle, whose role is to provide the energy supply 
necessary to bias the motion. Analytical expressions for the duration of the ATP 
cycle, for the Gibbs free energy and for the net displacement of the myosin head 
are obtained. Finally, by exploiting a method due to Sekimoto (1997, J. Phys. Soc. 
Jpn., 66, 1234), a formula is worked out for the amount of energy consumed during 
the ATP cycle. 
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1 Introduction 



The sliding process of a myosin head along the actin filament has been experimentally 
found by T. Yanagida and his co-workers to occur according to a random number of steps 
of constant size, mainly in a preferred direction . This observation is undoubtedly in strong 
support of the loose-coupling relation ( Oosawa . 2000[ ) existing between muscle contraction 
and ATP hydrolysis. The existence of such a loose coupling, in turn, assigns to the envi- 
ronm ental thermal energy an essential role. In a previous paper ((Buonocore and Ricciardi . 

it has been shown that a model based on a washboard-type potential can be de- 
signed in a way that, by a suitable tuning of the involved parameters, one can fit the avail- 
able data concerning the actomyosin interac ti on that is responsibl e for th e muscle contrac- 
tion m echanism (cf. iKitamura et al.l ( 19991 ): Ishii and Yanagida ( 2000f ): Kitamura et al. 
( 20011 ) and the references therein). 



The mechanism underlying actomyosin mechanical interaction as a result of the energy 
supply provided by ATP hydrolysis, and consequent transformation to ADP and Pi, has 
attracted a great interest not only within the realm of purely biological sciences but 
also among specia lists of stat i stical mechanics and soft condensed matte r. The variety of 
propo sed theories ([Magnasc o1 (|l99^ : lp"r"ost et al.1 (|l994l ): lAstumian1 (|l997l ): lMogimer et all 



(1998) and the references therein) share the aim of constructing models able to produce 
directional motion of a Brownian particle. 



In the present paper we construct a Brownian motor that embodies a fluctuating poten- 
tial ratchet within an idealized cycle of chemical reactions ensuing from ATP hydrolisis. 
The overall aim is to relate the net displacement of the myosin head to the amount of 
consumed energy in order to test whether the experimentally recorded displacements are 
compatible with the available total energy. While postponing to future papers data fit- 
tings and numerical evaluations, here we focus on working out analytical expressions for 
the duration of the ATP cycle, for the Gibbs free energy and for the net displacement of 
the myosin head. By exploiting the kinetic characterization o f the heat bath and the en- 
ergetics of thermal ratchet models along the lines indicated in Sekimotol (1997), a formula 
will also be obtained for the amount of energy consumed during the ATP cycle. 



2 The Model 



The aim of our model is to describe the dynamics of a myosin head during its inter- 
action with actin molecules. Here, the emphasis will be put on the construction of the 
model by means of rigorously quantitative arguments that will lead us to rigorous equa- 
tions and explicit formulas. The implementation of the obtained mathematical results for 
data fitting and prediction purposes will be the object of future endeavors. It should be 
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pointed out that, differentl y from previous papers (see Buonocore and Ricciardi ( 2003f ): 



Shimokawa and Sato! ( 200lh ) . our present approach focuses on energy features of a chemi- 



cally driven idealized gadget whose behavior is meant to mimic essential features of myosin 
head's dynamics and of related energy properties. 

Since the myosin motion occurs in a liquid environment, the effect of drag visco us friction 



must be included in the model. Moreover, on accounts of its size (see, for instance, Ishiwatari 



( 19971 )). a multiplicity of forces due to thermal agitation affect the dynamics of the myosin 



head, that in the sequel we shall often refer to as to "the particle" . Such forces, all ran- 
domly acting on the microscopic scale, are generally synthetized in a unique resulting ran- 
dom force. Hence, a random component will be included by us in the equation modeling 
the dynamics of the particle. Hereafter, we shall assume that the actomyosin interaction 
is of a conservative nature, and hence originating from a potential U(x,t), that we are 
viewing a space-time dependent. Denoting by x(t) the myosin head's position at time 
t under the customary point-size approximation, by m its mass and by p the existing 
viscous friction coefficient, the equation of motion is assumed to be 

dU(x,t) 



mx = -px 7^ \- ^2pk B T£(t) (1) 



where ks is Bolzmann constant, T is the absolute temperature and k B T (measured in pN • 
nm) is the thermal energy. The quantity 2pk B T originates from the fluctuation-dissipation 
theorem that also implies that is a zero- mean, delta-correlated Gaussian noise with 
unit intensity coefficient: 

<£(*)> = 0; (£(*)£(*')> = <*(*-*')■ (2) 



Eq. (1) alone does not suffice to describe the dynamics of the myosin head. Indeed, such 
dynamics is heavily biased by the chemical cycle elicited by ATP hydrolization, one of its 
effects being changes of potential U (x, t) in a fashion that is independent of the position 
of the myosin head. For our modeling purposes, we assume that ATP hydrolysis cycle is 
composed of three chemical states. These are depicted in Fig. 1, in which A, M, and Pi 
denote actin, myosin, and inorganic phosphate, respectively. 

It is worth pointing out explicitly that, as it is often customary when formulating a model 
of a complex phenomenon, the 3-state model proposed here is a drastic simplification of 
the underlying biochemical and biophysical mechanisms. Nevertheless, such simplification 
catches the essential feature of the dynamics of the real phenomenon. Indeed, it depicts the 
random "local" fluctuations as well as the abrupt changes corresponding to detachment 
and attachment of the myosin head. What is even more relevant, is that via our model we 
obtain exact formulas, which is an important propedeutic tool aiming at quantitative data 
fitting of the considered phenomenon. Furthermore, it should not pass unnoticed that the 
presently proposed model involves some twelve parameters (six transition rates, extreme 
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values, period and asymmetry parameter of the potential, viscous friction coefficient) so 
that it would be unrealistic at the present stage to increase further the complexity of the 
model if one wishes to make it suitable for fitting purposes. 



State 1 




Flashing 
ratchet 



State 2 


a23 


State 3 






n ; i 


■ 


a32 




M ■ ADP ■ Pi 




A • M ■ ADP • Pi 



How many steps? 



Fig. 1. The diagram of the chemical reaction route. 



The transition of the chemical state from state i = 1,2,3 to state j = 1,2,3 occurs 
according to a Poisson process with rates for i ^ j. The relations among transition 
rates Oy , equilibrium constants and equilibrium velocities k^ are indicated in Table 1. 

Table 1 

Equilibrium constants, velocity constants, and transition rates of the chemical reaction route. 



Equilibrium constant 


Velocity constant 


Transition rate (s 1 ) 


^12 = fa 2 /k 21 


k 12 (M-V 1 ) 
k 21 (M-V 1 ) 


a 12 = k 12 [ATP] 
a 2 i = k 2 i [A] 


K23 = &23A32 


k 23 (M-V 1 ) 
k 32 (s- 1 ) 


a 2 3 = hs [A] 
032 = &32 


K31 = hi/k\3 


hi (s- 1 ) 

ki 3 (M-V 1 ) 


031 = hi 

ai 3 = kis [ADP] [Pi] 



We assume that when the chemical state at time t is state i = 1,2, 3, then the potential 
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function U (x, t) is: 



U(x,t) = Ui(x) = < 



~(x-X a ), < x[mod A&] < X a 



(3) 



Ui 



, Xb — X a 



[x - X a ), X a < a; [mod A 6 ] < X b , 



where A;, is the length of an actin molecule and A a is chosen greater than A&/2 to break 
the symmetry of the potential function. We are thus invoking a periodic saw-tooth type 
potential to characterize each chemical state. When the particle moves from one potential 
well to a neighboring one, we conventionally say that the particle has moved one "step" . 
The net displacement of the particle at the end of an ATP cycle will thus be expressed as 
the difference between the number of steps moved in the positive direction and the number 
of steps moved in the negative direction, times the length Xb of one step. As sketched in 
Fig. 1, a relevant question relates to the number of steps during an ATP cycle. This figure 
also indicates that, at the end of the ATP cycle, one step in the prevailing direction is 
performed by the myosin head, which is a consequence of the allosteric conformational 
change following the release of Pi. Hereafter, A a , A& and all other involved lengths will 
be expressed in nanometers. The value Ui is the potential height for the chemical state 
i = 1,2,3. Since the binding strength between the myosin head and the actin filament 
is strong in states i = 1,3 and weak in state i = 2, we choose the potential maxima as 

u x > u 3 > u 2 . 



We now remark that it is conceivable to assume that the motion of the particle can be 
described under an overdamped approximation. Indeed, the viscous friction coefficient 
is given by p = 90 (pN • ns/nm), t he mass of the myosin head is m = 2.2 • 10~ 22 (Kg ) 



and Xh = 6 (nm) ( s ee, fo r instance, iBuonocore and Ricciardi l|200ah : | Ishiwatar 
Ishii and Yanagid"al (j2000j ). respectively). Hence, denoting by M = 100 (pN • nm) a sure 



upper bound to the potentials U±, U 2 , U$, one has 



m „, 2.2- 10" 1 -100 _ 4 
M = — -— ~ 10 -4 . 



p 2 A 



(9- 10-6) ; 



The non-dimensional quantity on the left-hand side b eing much smaller than 1, justifies 
the above-mentioned approximation ( Reimannl . 120 02). Since the motion is overdamped, 
the inertial component in Eq. (1) can be disregarded, so that the dynamics of the particle 
will henceforth be modeled by the following Langevin equation: 



x 




(4) 



where D = ksT/p is Einstein's diffusion coefficient. 
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Let now Z(t) = [X(t),I] be a random vector whose components X(t) and / have the 
following meaning: X(t) denotes the position of the myosin head at time t and J is a ran- 
dom variable, defined over the set {1, 2, 3}, denoting the chemical state of the actomyosin 
system. Then, 

fi(x, t) := ^-P {X(t) <x;I = t} (i = 1, 2, 3) 



is the joint probability density function (pdf ) of {X(t) G (x,x + dx) f] I — i}. Differently 
stated, fi(x, t)dx is the probability that the myosin head at time t is near x and that at the 
same time the actomyosin system is in chemical state i. The corresponding Fokker-Planck 
equation for fi(x, t) for the chemical state % — 1, 2, 3 is 



dt 



x,t) 



E a ij fi( X ' + E a kifk(x, t) 

i jjti I k^i 

\j = l ) k = l 



d_ 

dx 



l dUj 
p dx 



x)fi(x,t) 



+ D 



dx 2 



x,t).(5) 



For i = 1,2,3, 

+oo 

Pi{x,t) = Y fi( X + r) "b,t) 



is a periodic function of period A& that is a solution of Eq. (5) in which the space variable 
is made to vary in the interval (0, A&). Under such reduced dynamics, the existence o f 
a non trivial stationary regime is suggested by physical considerations ( Reimann . 12002^ . 
In this stationary regime, pf{x) = Pi(x) for states i = 1,2,3 are the joint pdf's of the 
position X of the myosin head and of the states I of the chemical system: 



Pi(x) = -^-P{X <x;I = i} 
dx 



1,2,3). 



(6) 



Such functions satisfy 







( 3 

E G * 



v. 



i=i 



3 d 
Pi{x) + a kiPk{x) + — 

fc=i 



-U-{x)pi{x) 
P 



+ Dp'!{x) 



1,2,3) 



or 



+ J2 a kiPk( x ) - S-i 



X 



1,2,3), 



(7) 



\i=i / 



k^i 
fc=l 
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where 



Si(x) = - -U[{x) + D^— Pi (x) (i = 1, 2, 3) 



(8) 



denote the probability currents. 



Note that adding term by term Eqs. (7) one obtains: 



S[(x) + S' 2 (x) + S' 3 (x) = 0, 



where prime means space derivative. Hence, 



S(x) := S^x) + S 2 (x) + S 3 (x) = S, 



(9) 



with S a constant denoting the total probability current of the myosin head at each 
x G (0, \b). Instead, 



is the pdf of the position of myosin head in (0, A^). 

We conclude this section by working out some results use of which will be made in the 
sequel. 

Proposition 2.1 Let g(x) be a periodic continuous function of period \ differentiable in 
(0, Aft) apart, at most, at an inner point \ a where, however, there exist bounded left and 
right derivatives. Then, 



Proof. Due to the assumed continuity and periodicity of g(x) one has g(0 + ) = g(0) = 

g(h) = g(K) and a(K ) = g(K )• Hence , 



p(x) := p x (x) + p 2 (x) + p 3 (x) 



(10) 




o 




J g\x) dx = J g'(x) dx + J g\x) dx = [g{\~) - g(0 + )} + [g(X;) - g(\+) 




= [g(K)-g^ + a)] + [g(K)-g(o + )} =o. 
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Corollary 2.1 Fori = 1,2,3 one has: 



a 6 

Jp' i (x)dx = (11) 
o 

At 

Js' i (x)dx = 0. (12) 



Proof. The proof follows immediately by nothing that the probability density func- 
tions (6) satisfy the assumptions of Proposition 2.1. Indeed, they are expressed as linear 
combinations of exponential functions in which the combination coefficients follow by im- 
posing that the Fokker-Plank equations are satisfied and that continuity conditions hold 
at the end points and A& and at the potential minimum X n . Simila r cond itions also hold 



for the probability currents. (See, for instance, lAstumian and Bierl f|l994r ) for a detailed 
argument). 



In conclusion, we point out that from Eq. (10) and from Eq. (11), that holds for i — 1,2, 3, 
one obtains: 

A„ 

p'(x) dx = (13) 

o 



3 State occupation probabilities and Gibbs Free energy 



The probability of the state i = 1,2,3 in the stationary regime is given by: 

A 6 

Pi = [ Pl {x)dx (14) 





which is independent of both x and t. Then, we have 

Proposition 3.1 The stationary state probabilities satisfy the following algebraic linear 
system: 

= -(a 12 + ai 3 )Pi + a 21 P 2 + a 31 P 3 

= a 12 Pi - (a 2 i + a 23 )P 2 + a 32 P 3 (15) 

= ai 3 Pi + a 23 P2 - (a 3 i + a 32 )P 3 
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Proof. System (15) immediately follows by integration of Eq. (7) for i = 1,2,3 with 
respect to x from to A& and by recalling Eq. (12). 



An alternative and more straightforward proof, that will also allow us to introduce some 
necessary notation and a useful result, will now be given. 

Proof. Let a = (a± 2 + a 2 \ + a 23 + 032 + 031 + a^)' 1 and let A be a positive quantity, 
with dimension of time such that A < a. Moreover, let P^ be the probabilities that 
the 3-state model undergoes a transition from state % to state j = 1,2,3) during A. 
The probability for the system to be in state i = 1,2,3 at time t satisfies the following 
equation: 

Pi(t + A) = P^Pu + P 2 (t)P 2i + Ps(t)P 3i . (16) 



In stationary condition, Eq. (16) yields: 



Pi — P\P\i + -^2-^21 + P 3 P 3 i — PiPii + X/ PjPji 

3=1 



or: 



o^iPu-m + y^PjP 



Jl- 



From the normalization conditions Pn + P^, + Pa = 1, for % = 1,2, 3, and from the remark 
that for i 7^ j one has 

Pij = A • Oij, (17) 



system (15) finally follows. 



The homogeneous linear system (15) has rank 2, so that it admits an infinity number of 
solutions. The states occupation probabilities are obtained via the normalization condition 
Pi + P 2 + P3 = 1. We further note that system (15) admits the invariant 

J c = a 12 Pi - a 21 P 2 = a 23 P 2 - a 32 P 3 = a 3 iP 3 - a 13 P 1 , (18) 
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representing the frequency of the ATP cycle in the 3-state model under the stationary 
regime. Its reciprocal Tq = 1/ Jc expresses the averaged time period of one ATP cycle. 
By multiplying both sides of Eq. (18) by A, one obtains another invariant expressed in 
terms of occupation and of state transition probabilities during time interval A: 

J C A = P 12 Pi - P 21 P 2 = P 23 P 2 - P 32 P 3 = P 31 P 3 - P 13 p. (19) 



We are now in the position to calculate Gibbs free energies AG i2 , AG 23 , AG 31 of all 
chemical reactions involved in the 3-state model. Indeed, assuming that the concentrations 
of the various reactants are proportional to the occupation probabilities of their respective 
states, one obtains: 



[M-ADP-Pi][A] P 2 [A] 
AG12 = —k B T In K 12 + k B T In — — — — — — = k B T m 



[A-M][ATP] X 12 Pi[ATP] 

= k B T\n^- (20) 
anPi 

AG23 = -* B T in A' 2 , + k B T ,n S^^gj = In ■ P » 



[M ■ ADP • Pi] [A] a K 23 P 2 [A] 

= k B T (21) 
d 23 P 2 

, , [A-MlfADPlfPil , PifADPlfPil 

AG 31 = —k B T In K 3 \ + k B T In = fc^T In ■ !! 



[A • M • ADP • Pi] P 3 A 31 

= A:bT1u^L (22) 
a 3 iP 3 

where use of the relations indicated in Table (1) has been made. 

Formulae (20), (21) and (22) explicitly show the dependence of Gibbs free energies on 
transition rates as well as on states occupation probabilities in the stationary regime. 

Proposition 3.2 Gibb's total free energy AG of an ATP cycle depends only on the tran- 
sition rates of the 3-states model. 

Proof. Adding sides by sides Eqs. (20) - (22) and recalling the relations in Table 1, 
one obtains 

AG = AG 12 + AG 23 + AG 31 = k B T In = k B T In ^i^H ( 23 ) 

Ai 2 A 23 A 3 i[ATPJ a 12 a 23 a 31 



10 



4 Step numbers and energy consumption 



The followin g proposition e xtends to the case of potential (3) the well-known property (see, 
for instance, Risken ( 1984 )) that the velocity and probability current are proportional. 



Proposition 4.1 The average speed of the myosin head in the 3-state model is given by 
v := E [x(t)\ = X b S (24) 



where S is the probability current defined in Eq. (9). 

Proof. Taking expectation on both sides of Eq. (4) and recalling the first of (2) we 
obtain: 



v = — E 
P 



dU(x,t) 
dx 



The expectation on the right-hand-side can be calculated by decomposition in terms of 
the three mutually exclusive and exhaustive chemical states. Hence: 



v = -~Y J P i E 

P i=l 



dU(x,t) 



dx 



-jlPi [ul(x)2& dx = --Y t fufaMx) dx .(25) 

P i=l n * P i=l n 



Recalling the definition (8) of the probability current when the system is in state % = 1, 2, 3, 
and making use of (9) and (10), from (25) we obtain: 



y 3 



\* 3 



8=1 



2 / [Si(x) + Dp'^x)] dx = JY, Si(x) dx + D J ^Tp'iix) dx 

i=l i=l ' 

S(x) dx + D J p'(x) dx . 
o o 

Use of Eqs. (9) and (13) complete the proof. 



From Proposition 4.1 we obtain the following 

Corollary 4.1 The net average number n of steps performed by the myosin head is given 
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by 



n 



S_ 



c 



(26) 



where S is the total probability current of the myosin head and Jc is the frequency of the 
ATP cycle in the 3-state model. 

Proof. The net displacement of the myosin head during one ATP cycle of duration Tc 
is 

v ■ T c = hS ■ T c = XbS ■ 



Dividing by the length of a step we are led to (26). 



Next, we shall calculate the mean energy consumed during one ATP cycle for the 3-state 
model, which depends on the potential differences among the states. When the particle 
occupies position x at t he time t when a transition takes place from state i to state j, the 



energy consumed is (see ISekimotol ljl997h ) 



H{x,t) = AU i: 



Uj{x) 



-UAx) 



(27) 



Let now r > denote an arbitrary instant. In (r, r + A) a state transition certainly 
occurs. The energy consumed during such transition is a function of (random) position 
X G (o, \ b ) of the particle at time r, of the (random) occupied state / G {1, 2, 3} in the 
3-state model and on the (random) direction of the transition I — > J, where J G {1, 2, 3}: 



H(X, t) = Hjj(X) = AU U (X). (28) 

Proposition 4.2 The average energy H during a transition occurring in (r, r + A) is 

given by: 



H:=E [H(X, r)] = P l P l2 {H l2 + H 21 ) + P 2 P 2 3(^23 + #32) + P 3 P 3 i(#3i + #13) + 

-J c A{H l3 + H 32 + H 21 ) (29) 



where 



H ij = E[H ij {X)}= H^x 



P 

1 7. 



A 6 



dx 



AU^^-dx 



(i,j = 1,2,3), (30) 
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the Pi 's are the stationary state probabilities and Pij 's are the probabilities of the transition 
i — > j in time interval A. 

Proof. A repeated use of averages of conditional expectations yields: 



3 3 

E [H(X, r)] =J2PiE [H(X, r) \ I = i] = y £P i E [H U (X)] 
i=i i=i 

=E^E p ^ I^M*) I j - = E ^ E ^ 

1=1 j = l j=l jr' = l 

3 3 

— E E PiP-ijHij- 

i=l j=l 

To complete the proof we notice that H u — for % — 1,2,3 and that, by virtue of the 
invariant (19), one has 

P 2 P 21 = P 1 P 12 - J C A; P 3 P 32 = P 2 P 23 - J C A; P 1 P i3 = P 3 P 3 i - J C A. 



We are now in the position to calculate the mean energy consumed during one ATP cycle. 
Proposition 4.3 The total mean energy Ti consumed during one ATP cycle is 

H = -r- [Pia 12 (H 12 + H 2l ) + P 2 a 23 {H 23 + H 32 ) + P 3 a 3l (H 3l + H 13 )] + 
■Jc 

-(H 13 + H 32 + H 21 ), (31) 

where = 1, 2, 3; i ^ j ) is the transition rate from state i to state j . 

Proof. The total consumed energy is expressed as a random sum of random variables, 
the number of terms to be added being the random number N of state transition that 
occur in the 3-state model during one ATP cycle. Furthermore, each such term is a random 
variable having the same distribution function as H(X, r). Hence, the mean total energy 
7i is obtained as the product of the mean value H of the random variable H(X, r) times 
the average number N of transition during one ATP cycle. Since the mean duration of 
one ATP cycle is Tc = 1/ Jc, we obtain: 



N:=E[N]=J2P t E[N\I = l ]j:P l j:^ = f E^E^ = f E^ = f 

i=l i=l j=l * i=l j=l i=l 
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Use of (17) and (29) completes the proof. 



5 Concluding remarks 

The necessity of relating the dynamics of the displacements of Myosin II heads during its 
interaction with actin filaments — as observed in the experimental work by T. Yanagida 
and his group — with the energy available during an ATP cycle has been the object of the 
present paper. To this purpose, we have designed a Brownian motor based on a fluctuating 
potential ratchet working as a gadget driven by a cycle of chemical reactions among 
hypothesized three different states of the actomyosin system. We have obtained rigorous 
analytical expressions for the duration of the ATP cycle, for the Gibbs free energies, for 
the net displacement of the myosin head and for the average energy consumed during an 
ATP cycle. Use of these theoretical results will be made in future endeavors in order to 
test in a quantitative fashion whether the existence of multiple steps is compatible with 
energy balance requirements. 
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